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Abstract 

Results in the quenched approximation for SU(3) breaking ratios of the heavy- 
light decay constants and the AF = 2 mixing matrix elements are reported. 
Using lattice simulations at Q/g 2 = 5.7, 5.85, 6.0, and 6.3, we directly compute 
the mixing matrix element M hi = (P h i\hj^(l - 75)^7^(1 - j5)l\Phl)- Extrap- 
olating to the physical B meson states, B° and B®, we obtain M^/My = 
1.76(10) in the continuum limit. The systematic error includes the errors 
within the quenched approximation but not the errors of quenching. We also 
obtain the ratio of decay constants, fbs/fbd = l-17(2)^t 6 2 . F° r the B param- 
eters we find Bb s (2GeV) = Bj,d(2GeV) = 1.02(13); we cannot resolve the 
SU(3) breaking effects in this case. 
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Using lattice methods, one can calculate the AF = 2 heavy-light mixing matrix element, 

MM = {P hl \h lp {l - j 5 )fh lp (l - 75) W • (1) 

As is well known, these matrix elements govern B° — B° and B® — B® oscillations In 
the above h and I denote heavy and light quark fields, Phi the corresponding pseudoscalar 
meson, and /1 is the energy scale appropriate to the calculation. Here we compute directly 
the SU(3) flavor breaking ratio [§], 

r sd = M bs (n)/M bd (n) . (2) 

Our central result is that r sd = 1.76(10)112 i n the quenched approximation, where the 
first error is statistical and the second systematic. The importance of this ratio is that, in 
conjunction with the eventual experimental measurement of B® — B® oscillations, it should 
allow the cleanest extraction of the crucial CKM parameter Vu- 

Since the CKM matrix elements are fundamental parameters of the Standard Model, it 
is clearly important to determine them precisely. V td is especially significant because low 
energy manifestations of CP violation, which enter through virtual t—t loops, invariably 
involve V td . At present, V td is deduced from B° — B° oscillations via the mixing parameter 
x bd = AM bd /T bd §. 

G 2 

Xbd = nd^m bd b{^)B bd {^)fl d M^rj Q cDS{x t )\V td \ 2 , (3) 

where m bd , r bd = T^ 1 , and f bd are the mass, life time, and decay constant of the B° me- 
son, and (AM) bd is the mass difference of the two mass eigenstates of the B° — B° system. 
x bd is the mixing parameter characterizing the oscillations and has been determined exper- 
imentally, x bd = 0.73(5) 0. B bd is the so called bag parameter, and b(fi) and S(x t ) are 
perturbatively calculated short distance quantities 0. To extract V td from Eq. ([$]) requires 
knowledge of two hadronic matrix elements, f bd and B bd . These are being calculated us- 
ing lattice and other methods. f bd may eventually be measured experimentally through, 
for example, the decay B — > tv t . However, B bd is a purely theoretical construct which 
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is inaccessible to experiment. Thus determination of Vtd from experiment will ultimately 
be limited by the precision of the nonperturbative quantity f bd B bd . These parameters are 
related to the matrix element Eq. ([[]) via 

MM = \furn 2 M B U: (4) 

and often one writes 6(/z)M^(/i) = M^, a renormalization group invariant (RGI) quantity. 

Making the replacement d — > s in Eq. (||]) and taking the ratio with Eq. (|3|), we arrive at 
an alternate way to extract Vtd, 

~\7/~~\2 = T *d (5) 

\Vts\ m bs r bd x bs 

Thus, in contrast to the above method for determining Vtd via use of Eq. (0), once the 
B° — B° oscillation parameter, x bs = ( A ^ ba is experimentally measured, we can use Eq. fl5|) 
to determine Vtd- The right hand side of Eq. (|5]) involves three SU(3) breaking ratios, only 
one of which, namely r s d, needs to be calculated non-perturbatively. The remaining two can 
be measured experimentally, at least in principle. Indeed, since the spectator approximation 
is expected to hold to a very high degree of accuracy J7| , it is also reasonable to expect that 
Tbs/^bd = 1 within a few percent. Of course, the measurement of x bs is very challenging. A 
variety of experimental efforts are underway at both e + e~ and hadronic machines towards 
that goal |§. Note also that V ts in Eq. (|5]) is related by three generation unitarity to V cb 
and is therefore already quite well determined, \V ts \ ~ \V cb \ = 0.041 ± 0.003 ± .002 0. The 
important distinction between using Eq. (|5]) instead of Eq. (|3|) is that the former requires 
only knowledge of corrections to SU(3) flavor symmetry while the latter requires the absolute 
value of the matrix element M b d- It is also important to realize that since r s d is a ratio of two 
very similar hadronic matrix elements, it is less susceptible to common systematic errors in 
lattice calculations, among which are scale dependence, matching of continuum and lattice 
operators, and heavy quark mass dependence. Indeed, the ratio r sd is, to an excellent 
approximation, RG invariant, even though the individual matrix elements M bs and M b d are 
scale dependent. 
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In passing, we recall that flavor symmetries have also played a crucial role in determining 
other CKM matrix elements. In particular, SU(3) flavor symmetry has been important in 
precisely determining V us = sin0 c . More recently, heavy quark symmetry (HQS) || has 
been used to improve systematically the determination of V c b. 

The lattice methodology for calculating these matrix elements (i.e. Eq. ([I])) is, by 
now, well known flO |. The amplitudes for Bq — Bq mixing, usually called "box" ampli- 



tudes, occur at second order in the weak interaction. After integrating out the W boson, 
the operator product expansion (OPE) allows one to write the corresponding amplitude 
as a short distance expansion. In this case there is only one operator in the expansion, 
Oll = &7/i(l — 7s)^>7/i(l — lh)d. Its (Wilson) coefficient, Cll(^), is calculated most easily 
in continuum perturbation theory. The matrix element of Oll must be calculated non- 
perturbatively on the lattice since it contains the long distance QCD information of the 
physical process in question. The product of the two yields the scale invariant amplitude, 
which is obtained by translating either result from one regularization scheme to the other. 
We accomplish this in the usual way by matching the lattice operator to the continuum 
operator in a particular scheme at some low energy scale. For convenience we choose the 
scale n = 2 GeV. Using the renormalization group equations, Cll(M w ) is then run down 
to this scale, which yields Cll{^)- 

For Wilson quarks the continuum-lattice matching for Oll has been carried out to one 
loop in perturbation theory ||il~|-|l3| . 



Ott = 4R h h (0% + yj^ ( Z + (a»)0% 



+ -^(20 ss + QOpp - llOw + UOaa + 20rr))) (6) 

where On corresponds to 7^(1—75) — * 1, 75, 7^, 757^, and in the expression for Oll- The 
Wilson quark action explicitly breaks chiral symmetry, so these new operators arise to cancel 
the chiral symmetry breaking terms in 1 ll- We use the naive dimensional regularization 
(NDR) scheme with "tadpole improvement," so Z + = (-50.841-4 In (a//) + 34.28) 0,|TT|jg 
where a is the lattice spacing. Z* = 9.6431 and depends only on the Wilson r parameter 



TABLE I. Summary of simulation parameters. Kh and ki are the heavy and light Wilson quark 
hopping parameters. 



Q/g 2 


conf. 


size 


flight 


Kheavy 


5.7 


83 


16 3 x 33 


0.160 0.164 0.166 


0.095 0.105 0.115 0.125 0.135 0.145 


5.85 


100 


20 3 x 61 


0.157 0.158 0.159 0.160 


0.092 0.107 0.122 0.130 0.138 0.143 


6.0 


60 


16 3 x 39 


0.152 0.154 0.155 


0.103 0.118 0.130 0.135 0.142 


6.0 


100 


24 3 x 39 


0.152 0.154 0.155 


0.103 0.118 0.130 0.135 0.142 


6.3 


100 


24 3 x 61 


0.148 0.149 0.150 0.1507 


0.100 0.110 0.125 0.133 0.140 


6.5 


40 


32 3 x 75 


0.146 0.147 0.148 0.1486 


0.100 0.110 0.120 0.132 0.137 0.142 



which we set to 1. The last term in Z + comes from mean field improved perturbation 
theory ||14|| , which removes tadpole terms. (Z* is an off-diagonal correction which does not 



have tadpole contributions at this order.) The scale at which the coupling g in Eq. ^| is to 
be evaluated is not fixed at one loop, however. It has been estimated for the decay constant 
using the methods of Ref. [0 as q* = 2.316/a jlB] . We use this scale to find the central 



values; the variation with two choices for the scale, 1/a and n/a, is used to determine the 
associated systematic error. The usual naive renormalization of the fermion fields, 4^^, is 



modified by the El-Khadra-Kronfeld-Mackenzie(EKM) norm [16] which is more suitable for 
the heavy quarks in our simulations. 

The Wilson quark action also introduces errors proportional to (powers of) the lattice 
spacing in observables. We attempt to remove these by simulating at several values of the 
coupling 6/g 2 (a) and extrapolating to a = 0. 

Table (|) summarizes the lattice data used in our analysis. For each K[ and Kh in Table | 
we calculate a quark propagator using a single point source at the center of the lattice and a 
point sink. These are contracted to obtain two and three point meson correlation functions 
which are fit simultaneously to obtain the matrix element Mhi- 

In Fig. p] we show sample results at 6/g 2 = 6.3 for Mm vs. k{ 1 for each value of Kh, where 
the quark mass in units of the lattice spacing is am q = (k^ 1 — k~ 1 )/2. Here, k c is the critical 



hopping parameter where the pion mass vanishes. Mm is extracted from the three point 
pseudoscalar correlation function, which is proportional to for large time separations of 
the four quark operator, and the two pseudoscalar meson interpolating operators. Results 
for the physical B and B s meson systems follow from a series of fits to the lattice data, 
which we use to extrapolate in the two parameters and K\. We use covariant fits and a 
jackknife procedure at each step to account for the correlations in the data. 

To begin the extrapolations, k c and k s (the strange quark hopping parameter) are de- 
termined from a fit to the squares of the pseudoscalar masses as a function of kJ 1 and /% 1 
(I and I' refer to non-degenerate light quarks ). We use the following fit form, which does 



not include the logarithmic terms relevant at very small quark mass [17 



m 



\, = C + Ci(K ; 1 + Kp 1 ) + C 2 K-i 1 K, l < 1 + C 3 (/^ 2 + K v 2 ). (7) 



A typical fit is shown in Fig. |2]. The values for k c and k s and x 2 /dof for each fit are 
summarized in Table [TT]. The curvature in Fig. [2] is small but certainly present: including only 
constant and linear terms in the fits generally yields poor x 2 values. The linear fits shown in 
Table |TT| were obtained by omitting the 3, 2, 0, 3, 8, and 6 heaviest points for 6/g 2 = 5.7 to 
6.5, respectively. The linear fits then had acceptable values of x 2 except at 6/g 2 = 6.0(24 3 ) 
and 6.3 where completely constrained fits were used. The values for k c obtained from the 
linear fits are in rough agreement with the quadratic fits; they are systematically low by one 
to two statistical standard deviations. Values for k c determined from the linear fits agree 



with earlier calculations flg-pO] at 6/g = 6.0 and |18j at 6/g = 5.85. At 6/V = 5.7, 



5.85, 6.0, and 6.3, k c is systematically higher by several statistical standard deviations than 
the values found by the MILC collaboration f2~I| , and the value at 6.0 in Ref. p2 |. In this 



study we use point sources on lattices with modest extent in the time direction. A detailed 



comparison with the data from Ref. [[H[ indicates that this is likely to be the main cause of 
the discrepancy. 

Since higher order chiral effects are completely different in the quenched and full theories, 
one might argue that the linear fits are preferable on physical grounds. For our central values, 



we stick with the quadratic fits, which describe our data well, but we take the difference 
arising from a switch to linear fits (as well as from the k c shift necessary to reproduce the 
Ref. [^TJ data) as an estimate of one source of systematic errors. 

Finding k s requires the scale a, which we set from af n , to determine the lattice value of 
the kaon mass amx is also tabulated in Table (0)). Our values for k s using linear fits 
agree to about one a with Ref. [|18|] at 6/g 2 = 5.85 and Refs. [H,|l^] at 6.0. Refs. fi8| , |l9H used 
m p to set the lattice spacing and amx, among others, to determine k s . Here we compare 
with values determined from amx- At 6/g 2 = 5.7, 5.85, 6.0, and 6.3, k s determined from 



arriK (with af w used to set the lattice spacing) agrees well with the results from Ref. |2T 



One might expect the values of k s , like k c , to disagree among the various calculations since 
they are determined from the same data. However, the added statistical uncertainty from 
the kaon mass is enough to mask the systematic error. We mention the above because the 
flavor breaking ratios given below are sensitive to the (relative) values of k c and k s . We also 
note that at 6/g 2 = 5.7 the choice of the coupling constant scale for Za, the lattice axial 
current renormalization which appears in the determination of f n , has a significant effect on 
the lattice spacing determination; Za differs by ~ 7% when the scale changes from 1/a to 
ir/a. 



Next, we linearly extrapolate Mm to Ki = k c and k s . The results for Mm at 6/g 2 = 6.3 (see 
Fig. [I]) show a smooth linear behavior. Similar results are obtained at the other couplings. 
Up to this step all of the covariant fits have acceptable values of x 2 , except the point at 
6/g 2 = 6.0 (24 3 ). Results at this point also showed significant variation with the form 
of the chiral extrapolation. The three point correlators here do not exhibit true plateaus 
but instead monotonically decrease with time, so there is undoubtedly contamination from 
excited states and additional uncertainty coming from the choice of fit range, which is 
necessarily small. Also, at 6/g 2 = 6.5 the data were too noisy to extract M^. Finally, the 
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TABLE II. Inverse lattice spacing and critical and strange hopping parameters. For each value 
of 6/g 2 , the two rows correspond to a determination of k c and re s by quadratic and linear fits, 
respectively, to the pseudoscalar spectrum. For the linear fits, the 3, 2, 0, 3, 8, and 6 heaviest 
points are omitted for 6/g 2 = 5.7 to 6.5, respectively. Each value of k c results in a corresponding 
lattice spacing from af n . % 2 /dof refers to the fit used to determine the quantity immediately to 
the left. An entry of "cf" means a completely constrained fit. 



Q/g 2 


o-^GeV) 


X 2 /dof 


K c 


X 2 /dof 




X 2 /dof 


5.7 


1.37(10) 


0.31/2 


0.16973(15) 


0.24/2 


0.1645(7) 


0.24/2 




1.35(9) 


0.31/2 


0.16953(9) 


0.43/1 


0.1640(8) 


0.43/1 


5.85 


1.65(13) 


0.01/1 


0.16170(8) 


0.33/3 


0.1576(6) 


0.33/3 




1.64(13) 


0.01/1 


0.16157(5) 


1.30/4 


0.1575(9) 


1.30/4 


6.0(16 3 ) 


2.03(17) 


3.16/1 


0.15725(23) 


0.62/1 


0.1544(5) 


0.62/1 




2.01(16) 


3.16/1 


0.15715(6) 


1.08/3 


0.1545(4) 


1.08/3 


6.0(24 3 ) 


2.08(13) 


0.67/1 


0.15714(4) 


2.6/1 


0.1544(4) 


4.0/1 




2.17(15) 


0.67/1 


0.15739(4) 


cf 


0.1548(4) 


cf 


6.3 


3.09(21) 


0.81/2 


0.15199(4) 


9.5/6 


0.1502(2) 


9.5/6 




3.10(21) 


0.81/2 


0.15191(4) 


cf 


0.1503(2) 


cf 


6.5 


4.29(50) 


0.38/2 


0.14993(18) 


1.04/5 


0.1486(3) 


1.04/5 




4.22(49) 


0.38/2 


0.14972(1) 


3.04/1 


0.1487(3) 


3.04/1 



point at 6/g 2 = 5.85 is somewhat problematic. The statistical errors are large, so this point 
does not have a large impact on the continuum extrapolation. The difficulty arises in the 
three point correlators which show plateaus with a somewhat large oscillation. Three of the 
four light re's happen to be below k s while the fourth is just above. Thus, our light re's are 
closely spaced. The above considerations lead to a relatively inaccurate determination of 
the slope of M hl vs. re; -1 , which essentially determines r sd . 

The heaviest mass points in our calculation suffer from heavy quark systematic errors; the 
lattices are too coarse to resolve objects with mass greater than the inverse lattice spacing. 
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The biggest correction of these errors comes from using the EKM norm mentioned above. 



An additional correction can be made by using the so-called kinetic mass JT6[ in place of 
the meson (pole) mass in the heavy mass extrapolations described below. As in Ref. |23| . 
we use the tadpole improved tree level definition of the kinetic mass, m k P m = mp + rfi2 — m. 
m 2 and m are the tadpole improved heavy quark kinetic and pole masses, respectively. This 
definition is motivated by a non-relativistic expansion of the heavy-light meson mass and 
reduces to the usual meson pole mass in the limit where the heavy quark becomes light. 
For the heaviest masses, the kinetic pseudoscalar mass is almost double the pole mass. This 
correction is also used in Refs. [pl.EHl. 



We fit M hc to the HQET form 

M hc = C^ilTLp + C + Ci — . (8) 

m P 

Here mp is any definition of the heavy-light pseudoscalar mass. The resulting fit is evaluated 
at the experimentally known B° meson mass to determine the physical value of the matrix 
element. For the heavy-strange case we first extrapolate Mm to k s instead of k c . The form 
in Eq. |^ follows from the HQET results for the decay constants []T],|2|,|S| and the B parameters: 

fpVmp = d + dx/mp + 0(mp 2 ), (9) 

Bp = b + b 1 /m P + O{m P 2 ). (10) 

Our data are consistent with these forms. We note that for each value of6/g 2 all of the data 
points are covariantly fit to the above form, and each fit has a good confidence level except 
the one at 6.0 (24 3 ). An example is shown in Fig. [3[ It is noteworthy that the data fit the 
form in Eq. [8] over such a large range (this is true for all of the couplings we studied) . At 
each coupling the heaviest (kinetic) mass is close to the physical B mass. 

Fig. (|]) shows r s d = M^/Mm as a function of a. The ratio is greater than unity for each 
value of Q/g 2 . Using Eq. [], the renormalization scale is set to // = 2 GeV and the coupling 



is evaluated at q*. r s d is also tabulated in Table [HI . 



TABLE III. Summary of results for r s d- The last two rows refer to constant and linear con- 
tinuum extrapolations, respectively. The errors shown in parentheses are statistical. Column 1 
gives central values; columns 2-8 represent systematic differences in r S( i relative to column 1 and 
are used to estimate the corresponding error (see text). 



6/a 2 


1 


2 


3 


4 


5 


6 


7 


8 


5.7 


1.65(15) 


1.64 (15) 


1.50 (15) 


1.52 (25) 


1.52 (25) 


1.68 (15) 


1.59 (14) 


1.76 (20) 


5.85 


1.78(25) 


1.75 (24) 


1.44 (30) 


2.87 (2.20) 


2.07 (45) 


1.77 (27) 


1.68 (21) 


2.21 (27) 


6.0(16 3 ) 


1.80(23) 


1.79 (23) 


1.64 (25) 


1.70 (23) 


1.76 (22) 


1.74 (19) 


1.66 (19) 


1.79 (23) 


6.3 


1.96(23) 


1.94 (23) 


1.95 (33) 


1.85 (34) 


1.97 (26) 


1.83 (19) 


1.68 (16) 


1.99 (30) 


oc 


1.76(10) 


1.75 (10) 


1.58 (8) 


1.67 (15) 


1.77 (13) 


1.74 (9) 


1.64 (8) 


1.90 (12) 


oc 


2.18(39) 


2.21 (42) 


2.03 (42) 


2.10 (56) 


2.29 (48) 


1.94 (34) 


1.76 (29) 


2.09 (49) 



As mentioned earlier, we expect the Wilson quark action to introduce discretization 
errors of order a in all observables. However, for the ratio of two similar quantities, we also 
expect a significant cancellation of these errors. A constant fit gives M^/Afy = 1.76(10) 
while a linear extrapolation in a gives 2.18(39) (column 1 in Table p|) . The above fits have 
small x 2 values due to the large statistical errors, and we cannot rule out one fit in favor 
of another based on \ 2 ■ The measured slope for the linear fit differs from zero by ^ la. 
The linear trend may easily disappear with a one standard deviation change in either of the 
two end points, so we use the constant fit as our central value and the linear result as an 
estimate of the systematic error in the continuum extrapolation. 

Next we estimate other systematic uncertainties in our analysis. The details are given in 



Table [ITJ. Columns 2-8 refer to separate analyses where one parameter was changed from its 
reference value used to obtain column 1. The difference in the new extrapolated value is then 
taken as an estimate of the systematic error in r S( ^. In the following we list the uncertainties 
(numbers in parentheses refer to the corresponding column in Table |TJ). (2.) Changing 
the coupling constant scale to a" 1 yields r s ^ = 1.75(10) and 2.21(42) for constant and linear 
continuum extrapolations, respectively. (3.) Using quadratic chiral extrapolations for the 
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matrix elements yields 1.58(8) and 2.03(42). The fits used for the chiral extrapolations were 
completely constrained except at Q/g 2 = 6.3 and 5.85, each of which had one degree of 
freedom. (4.) Using the uncorrected pole mass yields r sc t = 1.67(15) and 2.10(56). For this 
last case we note that x 2 values were uniformly poor. The resulting fits underestimated the 
data at the heavy masses; the lighter masses, which had smaller statistical errors, dominated 
the fits. (5.) Constrained fits using only the heaviest masses yield 1.77(13) and 2.29(48), 
which we use to estimate the systematic error from including heavy masses that may be too 
light. (6.) As previously noted, we expect r S d to be sensitive to small relative shifts between 
k c and k s . Using the linear extrapolations for k c and k s , we find 1.74(9) and 1.94(34). The 
main effect is to lower the value of r s d at Q/g 2 = 6.3, which primarily affects the linear 
continuum extrapolation. (7.) As mentioned earlier, our values of k c may be systematically 
high. Shifting k c by -0.0003 at each coupling gives r S( i = 1.64(8) and 1.76(29). The shift was 



estimated from the difference of our k c values with those of Ref. where Gaussian sources 
and longer lattices in the time direction were used. While the absolute shift is numerically 
small, it amounts to several statistical standard deviations and is thus not accounted for in 
the jackknife analysis. 

(8.) Finally, we estimate the systematic error resulting from changing the fit range 
of the three point correlation functions. Thus far acceptable values of \ 2 were obtained 
using covariant fits to the three point correlation functions. However, with point sources 
(and lattices with modest extent in the time direction) the correlators do not exhibit long 
plateaus, so the allowed fit range is necessarily small. The fit ranges were shifted up or down 
by one or two time slices at each 6/g 2 which generally resulted in worse \ 2 fits- The only 
appreciable variations were at Q/g 2 = 5.7 and 5.85. 

Columns 2,3,4,6, and 7 give lower results than the central value (column 1) for the 
preferred (constant) continuum extrapolation. Adding the differences linearly gives a sys- 
tematic error estimate of —0.42. Combining the positive differences in columns 5 and 8 
with that from the linear continuum extrapolation gives +0.57. The final result is then 
r sd = 1.76(10)jj£. 
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After our initial determination another group calculated r s d in the static approxima- 
tion at Q/g 2 = 6.0. They find 1.35(5) pB| . When extrapolated to the static limit, m~p l = 0, 
our data at Q/g 2 = 6.0 yield r s d = 1.39(30) which agrees well with the above. In addition, 



our data at each value of Q/g 2 indicates that r s d is a smoothly increasing function of 



in 



p ) 



thus the static result may be a lower bound for r S( j. 

The extraction of the individual values of Mm and M bs is clearly expected to have 
larger errors. Thus, conventionally H^Jl|,§] these matrix elements are given in terms of 
the corresponding B parameter, which is better behaved. Carrying out a similar continuum 
extrapolation as above for Bbd{fJ>), we find a constant fit yields 5^(2 GeV) = 0.97(3) while 
linear extrapolation in a gives 1.02(13). We cannot, however, distinguish Bb s (2GeV) from 
Bbd(2 GeV) since our data for B^i vs. kJ 1 are fit equally well to constant or linear fits. This 
was not true for Mjy, as is evident from Fig. [|. Using linear extrapolations in both a and 
we quote 5^(2 GeV) = _Bb s (2GeV) = 1.02(13), where the error is purely statistical 
(systematic errors are small in comparison). 

We recall that until now lattice results for the SU(3) breaking ratio r s d have been 
obtained by using Eqs. (Q) and and the lattice determinations of fbd(s) an d B hd ( s )- A 
simultaneous fit of the pseudoscalar and axial vector correlators yields the decay constant 
fhi- Using Eq. § plus corrections up to 0(m^ 2 ), we find for the ratio of B meson decay 
constants, f bs /f bd = 1.17(2)(+2)(+5)(+l)(-2)(±4) (see Fig. | and Table 0). This result 



is for a constant continuum extrapolation which is reasonable for the data shown in Fig. |5|. 
The uncertainties are statistical, and the following systematic differences from: (2.) using 
the pole mass instead of the kinetic mass, (3.) using quadratic chiral extrapolations for the 
heavy-light mesons instead of the linear ones, (4.) replacing quadratic chiral extrapolations 
with linear extrapolations in the determination of k c and k s , and (5.) a shift in k c as before. 
Again, the numbers in parentheses correspond to the columns in Table 0. Finally, we 
consider the overall uncertainty in the slope of the ratio versus lattice spacing. A linear 
continuum extrapolation using all the data has a negative slope; while omitting the point 
at Q/g 2 = 5.7 yields a positive slope. Also, the results in columns 4 and 5 have positive 

12 



TABLE IV. Summary of results for the ratio fbs/fbd- The last two rows refer to constant 
and linear continuum extrapolations, respectively. The errors shown in parentheses are statistical. 
Columns 2-5 represent systematic differences relative to column 1 (see text). Where there is no 
entry, a reasonable fit was not found. 



r 1 2 


1 

1 


2 


3 


4 


5 


^ 7 




1 Ifi^i (Q7) 

± . lUu \ iy t J 


1 (p.(\\ 

i.iui luu j 


1 181 (10) 

l.lul l OU I 


1 14D (OK) 


5.85 


1.190 (80) 




1.493 (249) 


1.184 (79) 


1.287 (95) 


6.0(16 3 ) 


1.187 (36) 


1.180 (92) 


1.204 (61) 


1.182 (34) 


1.172 (35) 


6.3 


1.159 (32) 


1.180 (93) 


1.266 (74) 


1.177 (33) 


1.149 (31) 


6.5 


1.190 (66) 


1.289 (192) 


1.408 (243) 


1.145 (40) 


1.100 (46) 


oc 


1.167 (17) 


1.184 (52) 


1.214 (37) 


1.174 (17) 


1.148 (16) 


oo 


1.183 (55) 


1.241 (169) 


1.375 (138) 


1.152 (45) 


1.133 (44) 



slopes. Thus we include a symmetric error of ±0.04 due to the continuum extrapolation. 
F° r hs/hd) there were no significant differences due to changing the fit range by one or two 
units. Adding the above systematic errors, we find hs/fbd = 1.17(2)1^. We have omitted 
the larger volume at Q/g 2 = 6.0 from the analysis for reasons similar to those described 
above. In addition, the data set at 6/g 2 = 5.85 causes the same difficulties as before. For 
the ratio of decay constants we are able to get a statistically significant result at 6/g 2 = 6.5 
which is included in the above analysis. 

Our result for fbs/fbd is consistent with previous estimates |l|,|2l j23|j28[| . Note that while the 
decay constant using Wilson quarks has a perturbative correction (which does not depend 
on the scale fi), it cancels in the ratio (up to small quark mass corrections). As indicated 
above, the ratio of B parameters is consistent with unity, and the ratio of masses is 1.017 
0, so the old method leads to r s d ~ 1.42(5)ii5 which is compatible with, though somewhat 
lower than, 1.76(10)^42 from our direct method. As we have emphasized, the direct method 
has many desirable features which may allow future lattice computations to significantly 
improve the precision of this method for the determination of the ratio r s ^. 
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6.6 6.7 6.8 

FIG. 1. The four quark matrix element Mm at 6/g 2 = 6.3 as a function of the inverse light 
quark hopping parameter. Results for the other values of 6/g 2 are similar. The solid lines are 
covariant linear fits to the data. 
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FIG. 2. The pseudoscalar mass squared as a function of the non-degenerate light quark 
hopping parameters. 6/g 2 = 6.3. Results for the other values of 6/g 2 are similar. The solid lines 
are from a covariant fit to the form in Eq. [?]. 
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0.2 0.4 0.6 0.8 

nCCGeV" 1 ) 

FIG. 3. The four quark mixing matrix element for down(octagons) and strange(squares) pseu- 
doscalar mesons as a function of the inverse heavy-down(strange) meson mass at Q/g 2 = 6.0(16 3 ). 
The solid lines are covariant fits using Eq. [8| to all of the data points. We find very similar results 
at the other values of 6/g 2 . 
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FIG. 4. The SU(3) flavor breaking ratio r S( i = M^/M^ vs. the lattice spacing a. /i = 2 GeV 
and the coupling in Eq. || has been evaluated at q*. The lines denote constant and linear fits to 
the data, fancy squares the corresponding continuum extrapolations. 
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FIG. 5. The SU(3) flavor breaking ratio of decay constants fbs/ fbd vs - the lattice spacing a. 
Results are plotted for linear chiral extrapolations. The line denotes a constant fit, and the burst 
is the corresponding continuum extrapolation. 
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